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We consider the impact of second order corrections to the geodesic equation governing gravi¬ 
tational lensing. We start from the full second order metric, including scalar, vector and tensor 
perturbations, and retain all relevant contributions to the cosmic shear corrections that are second 
order in the gravitational potential. The relevant terms are: the nonlinear evolution of the scalar 
gravitational potential, the Born correction, and lens-lens coupling. No other second order terms 
contribute appreciably to the lensing signal. Since ray-tracing algorithims currently include these 
three effects, this derivation serves as rigorous justification for the numerical predictions. 


I. INTRODUCTION 

Gravitational lensing of background galaxies by large scale structure offers an excellent way to study the distribution 
of matter in the universe Emil Measurements of the cosmic shear are already enabling us to constrain the dark 
matter abundance @,|dl. In the future, large surveys may well uncover properties of dark energy, such as its abundance 
and equation of stateQ, ism and of neutrinos llll . This program will be successful only if we can make very 
accurate predictions about what theories predict El- 

Predicting the cosmic shear signal is more difficult than making predictions for the cosmic microwave background 
but much more straightforward than for the galaxy distribution. The former rely predominantly on linear perturbation 
theory so are extremely robust. The latter require not only nonlinear evolution but also an understanding of how the 
galaxy distribution is related to the mass. For lensing, we need to account for nonlinearities, i.e. move beyond linear 
perturbation theory, but, since the light trajectories depend only on the mass, we do not need to worry about how 
and where galaxies form 1 . 

The obvious way to treat nonlinearities is to perform numerical simulations of structure formation. Since dark matter 
is much more abundant than baryonic gas, only N-Body simulations are needed, not 2 the much more computational 
costly hydrodynamical simulations. A number of groups have run such simulations and extracted predictions by 
ray-tracing photon paths through the simulated universe El El In principle, this proceedure accounts for some 
second order effects, such as the second order perturbations to the potential and the fact that the path of the photons 
is not simply a straight line. However, the ray-tracing algorithims use the first-order geodesic equation to find the 
distortions in each slice of the universe. Therefore, again in principle, the algorithms may be missing relevant terms. 
A fully consistent second order calculation is warranted. 

In this paper, we start from the second order metric, proceed to obtain the second order geodesic equation, and 
then compute the cosmic shear consistently retaining all second order terms. The rest of the paper is organized as 
follows: §11 reviews the first order calculation, and §111 presents an order of magnitude estimate that will prove useful 
when wading through a host of second order terms. §IV presents the full second order geodesic equation, borrowing 
heavily from previous work. In it, we winnow out the terms that contribute negligibly, keeping only those terms that 
are responsible for observable deflections. 


II. FIRST ORDER SHEAR 

Consider the first order deflection of light due to inhomogeneities along the line of sight. We work in a flat universe; 
then the geodesic equation may be integrated to give 
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1 We refer here to cosmic shear. Other lensing phenomena, such as galaxy-galaxy lensing, do involve the interplay of light and mass. 

2 This ceases to be true on very small scales where the effects of baryons can no longer be neglected El 13- 
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Here x± is the perpendicular deflection of a ray starting at comoving distance away Xs and traveling towards us at 
X = 0 detected at angular position 9 with respect to a fixed z-axis. The first superscript on x±_ denotes the order of 
the perturbation, the second is a space-time index; we are interested only in the two components perpendicular to 
the line of sight. The first order distortion is 


/i 1}< 


M = -r SV 0)C 


MP 


( 2 ) 


Here T^ is the first order Christoffel symbol and the zero-order space-time vector p(°) a = (1, — e 1 ) ~ (1,— 9 1 ,— 1), 
where e l is the spatial direction vector. The square of e is unity, and 9 1 is the transverse displacement from the 
z-axis. We work in the small angle approximation throughout, so 9 is assumed small, and the z-component of e is 
approximately unity. 

In the contraction of the Christoffel symbol with the two factors of p, the only terms which contribute are those 
with a, (3 = 0 or 3. Using standard results 0, we can write the contraction as 


p( 1 ) I r ,(o)“ rj (o)/3 
1 a0 P P 
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where q^ is the first-order gravitational potential evaluated at the unperturbed 3 position of the light: £ (0) ( 9 , x) = 
X[M]- So, 


x { l ]i = -2 
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The shear matrix is defined as the derivative of this perpendicular deflection with respect to observed angle 9: 

1 dx\ 
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To lowest order of course x^ 1 = Xsd z , so the zeroth order term in the shear matrix is simply the identity. From 
Eq. Q, we see that the first order term in the shear matrix is 
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Again the second equality here follows because the potential is evaluated along the unperturbed path of the light. 


( 6 ) 


III. ORDER OF MAGNITUDE ESTIMATE 

In order to wade through the second order terms and extract the most relevant ones, we need a way of doing order 
of magnitude estimates. From Eq. ©, the shear induced by a perturbation with wavelength A is of order (rjj/ A) 2 </> 
since cosmological distances x are of order the Hubble radius, rn = 3000h -1 Mpc. The mean of the shear of course 
vanishes, but its rms we would expect to be of order (r£r/A max ) 2 </> rms , where A max is the wavelength near which 
perturbations contribute the most to the deflections. In currently popular models, A max ~ 30h 1 Mpc. The rms 
amplitude of the gravitational potential is of order 10~ 5 , so one would naively expect the rms shear to be about ten 
percent. The rms is actually smaller than this because not all Fourier modes contribute to the variance. Perturbations 
which vary rapidly in the z-direction lead to little total distortion since regions of positive and negative overdensity 
along the line of sight cancel each other. Thus, in the 3D Fourier space only modes with very small (< r^ 1 ) 
contribute. The fraction of modes which satisfy this constraint is (A /rn) for a given A [I, j]; the variance is therefore 
smaller than the naive estimate by this fraction, and the rms by its square root. The rms amplitude is therefore of 
order (rjy/A max ) 3/,2 (/> rms , less than a percent. 


Evaluating it at the perturbed position leads to a second order term, and for now we are considering only first order terms. 
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There are two lessons we learn from this order of magnitude estimate which will be useful when we evaluate second 
order terms. First, we need consider only those terms in /O which vary little along the line of sight. Terms with partial 
derivatives with respect to z = x 3 can be neglected. Second, the shear is much larger than the depth of the typical 
potential well, </> rms - It is the rapid changes in the potential which lead to large deflection; i.e., (d(f>/dx l )r h </> (as 
long the derivatives are in one of the transverse directions, i = 1,2). When considering second order terms, therefore, 
we will be most impressed by those with the most derivatives. 

Armed with this information, let’s return to Eq. m for /[' *’ and consider the correction incurred when we account 
for the fact that the photon does not travel along a straight line so the argument of (j> might reasonably be taken as 
x ^ + x W. Expanding about the zero order path then leads to one possible second order term 


d 2 (j) 

dxidxj 


x Wj 
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This term is known as the Born correction\vA ITsj . We will soon see that it emerges as one of many second order 
terms. To derive it here, we have cheated since the argument of <fi in Eq. m is the zero order path. We introduce 
this term now only because we know it will show up in the second order zoo, and we want to estimate its order of 
magnitude. From Eq. 0, is of order so the Born correction is of order (d 2 (f>)r^ I d(j). Thus, we need keep 

second order terms in only if they are of order r\d V- 


IV. SECOND ORDER GEODESIC CORRECTIONS 

The perpendicular deflection to second order is given by the analogue of Eq. 0 , with the superscripts changed 
from O to and [19] 

f(2)i = _ r ( 0 )i a/ 3 p (i)ap(i )/3 _ 2r ( 1 )i a/ 3 p ( 0 ) “p (1)/3 - 2d CT r ( 0 )i a/ 3 x ( 1 ) > ( 0 ) V 1)/3 

- - ^dvdrT^apx^x^p^p^ 0 - T™* af3 p {0)a p i0)t3 . (8) 

The zero order and first order Christoffel symbols are well-known mu, while the second order has been 
computed by Bartolo et al. (2l| . The zero order direction vector is given after Eq. m, while the time and space 
components of the first order direction vector are 


p ( 1 ) 0 (x) = p (1)0 (x = 0 ) - f / ( 1 ) 0 (xW, 

Jo 

p {i) \x) = - r 

Jo 


(9) 


with /W a = _r( 1 )“ (L4J ,p(°)A l p( 0 ) 1 '. To compute the first and second order Christoffel symbols, we need to specify the 
metric, thereby choosing a gauge. Following the formalism of Jsfi}, we set 
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goi = a 2 
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1 + 2 ^ 

1 , ( 2 ) . 1 ( 2 ) 

-diU ( 1 + -u\ 


^1 - 2^ (1) 


&U + \ D ax (2) + \ (dixf ] + 0jXi 2) + Xif 


( 10 ) 


( 2 ) ( 2 ) ( 2 ) 

Scalar perturbations are described by (j >, i/j, x an( I vector perturbations by lo] and x\ i an d tensors by Xij ■ 
Finally 


Dij = didj - ^V 2 . (11) 

Note that we assume that there are no first order vector or tensor perturbations. Even without these, scalar, vectors, 

( 2 ) ( 2 ) ( 2 ) 

and tensors mix at second order, so x\j an< I an d X) are generally non-zero, quadratic in first order scalar 
perturbations. 

We now work through the terms in Eq. © explicitly, retaining only those of order rj^d 3 ^ 2 or higher. 










First Term: — r( 0 ) l aj gp( 1 )“p( 1 )^ 

The zero order Christoffel symbol is proportional to the Hubble rate and is non-zero only when one of the lower 
indices is equal to i and the other equal to zero. Thus this first term reduces to 

- r (0)i a/3 p (1) > (1)/3 = -2r^V 1) V 1)i - (12) 

Forgetting about the boundary term in the first of Eq. ©, and remembering that \ ~ tr, we see that pH) 0 ~ 
Cst/^ 0 andp) 1 ) 1 ~ So this first term in /) 2 ) 1 is of order tr f ^ 1 ■ We showed in §11 that /H)* = 2(f) ^ 

so this first term is of order rird</>/H) 0 ; it contributes appreciably only if /H)° is of order rnd 2 (f) or higher. But 

/ (1) 0 = -r^°p (0) > (0)/3 . (13) 

Recall that p^ a is negligible unless a = 0 or 3. If both a and (3 are zero, then the Christoffel symbol is 
<j> ~ rjf(f> <C rnd 2 ^. If both indices are equal to 3, the Christoffel symbol again is of order 4>/th, he. negligible. 
If one of the indices is zero and the other equal to three, the Christoffel symbol is equal to 83 <f>, which is negligible 
for all perturbations of interest, i.e. for all perturbations which vary little along the line of sight. Thus this first 
term does not contribute. 

Second Term: — 2T^p^ a p ( ' 1 ' )/i 

Since only the a = 0 or 3 components of p)°) are non-negligible, 

- 2r2V 0) > (1)/3 -► - 2 r$V 1)/5 + 2 rgV 1)/3 ( 14 ) 

Recall from the previous paragraph that 

pH) 0 rv, (f) 

p {1)l ~ r H di<j) (15) 

So the biggest contribution will come from (3 = j, with j = 1 or 2, one of the transverse directions. The 
contribution is of order rud(f>T. But T^ 1 = ~ r^ 1 ^, so the first term on the right in Eq. is of order 

8 (f ) 2 and can be neglected. The second term in Eq. 11411 . 2 ry l ) , pHb , is of order r^dipT^j 1 ; recall 0 that this 
Christoffel symbol sets rotating pairs of its indices equal to each other, with the other index applying to the 
derivative of the potential. We certainly then do not want the index 3 to be the derivative since 83 is very small. 
But both i and j are transverse indices so cannot be equal to 3: all terms here vanish. 

Third Term: — 29 CT r( 0 ) i a p:r) 1 ) CT p) 0 ) Q p) 1 ) /3 

After invoking some of the approximations used in the previous paragraphs, we get 

- 2 d a T^ i a 0 x Wa p Wa p {1)/3 -> -28 a T^ 0j x Wa p Wj . (16) 

The derivative here acts only on the Christoffel symbol, which depends only on time. Thus a = 0, and cT ~ rj£. 
The first order perturbation in a,’H)° is of order rR(f >, so this term is of order defy 2 , far too small to contribute. 

Fourth Term: — d o Tn)*„p£( 1 ) cr p( 0 )“p(°)' S 

This is the only term in Eq. © quadratic in the first order variables (i.e., the only one of the first five terms) 
which contributes. It gives the Born correction and lens-lens coupling, a very nice result emerging naturally 
from this second order formalism. To see this, note that the only terms in the implicit sum which contribute 
are those with a = /3 = 0ora = /3 = 3. Both of these contribute identically, so this term is 

- 2 = - 4 aS? £ d *'<* - *'>!?■ < 17 > 


Fifth Term: -{l/2)8 a d T T^ i a px^> a x^ T p^ a p^ 

This term vanishes since the zero order direction vectors are non-negligible only if the indices a and (3 are set 
to 0 or 3, and r(°)' 00 = r(°) i 03 = r(°) i 33 = 0. 
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• Sixth Term: — T^ 1 a pp(°) a p(°)d 

The only remaining term is the one explicitly second order in T. Before delving into this term, we can localize 
it further by noticing that the indices a and (3 must be either 0 or 3. Thus, this last term is 


r (2 )\)o+ ^33-2^03 


(18) 


To compute T^ 2 ), we lift the results from the Appendix of Ref. [2]}. There, they compute the Christoffel symbols 
to second order. There are a number of simplifications we can make. First, we work in a gauge (the so-called 
Poisson gauge) in which the scalars u/ 2 ^ and x*' 2 '* as well as the vector x^ are set to zero, so vector perturbations 
are described only by u/ 2 ) 1 , henceforth simply called u>. Similarly, we assume there are no tensor perturbations 
at first order, and the second order perturbations are described solely by the traceless and divergenceless tensor 
Xij [spatial indices are raised and lowered with the Euclidean metric, so we take no care in distinguishing upper 
from lower indices]. Within this framework, 




__ - c,i 

00 T 1 33 03 — 2 U 


^(2) + ^(2) + 40 (1) a z </> (1) + - [uj 11 - d 3 u > 1 + diU 3 ] + - [2d 3 Xi3 - diX 33 - 2x' 3 ] ■ 

(19) 


Here prime denotes differentiation with respect to conformal time. 


To weigh the relative importance of each of these terms, it is important to remember that the only terms that 
contribute are those which are of order r 2 H d 3 cj) 2 . This results in both vector and tensor modes being irrelevant. 
Let’s work through the vector case explicitly; the tensor case is similar. The vector modes have been worked 
out in Ref. fH}. Eqs. (17-19) there give explicit expressions for the second order vector perturbations (their 
V = —u>/2). By counting powers of k —> <9, we see that to 1 ~ rndcfr. Even if we forget about the constraint that 
d 3 must be small, this means that the second line is at most of order r#<9 2 0 2 , too small by a factor of order 
th/X ma x to contribute. A similar argument applied to the expression for tensor modes in Ref. {22j leads to a 
similar conclusion. Thus, the only relevant explicit second order terms are from scalar perturbations, those on 
the first line of Eq. m- Using the same order of magnitude estimate, we can immediately eliminate the second 
scalar term on the first line, since it is of order <9</> 2 , clearly too small to be relevant. 

The only relevant explicit second order terms therefore are 


-p(2)z 
1 00 
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The final limit comes from recognizing that anisotropic stress, inducing a iji — </> ^ 0 term, is a post-Newtonian 
second-order contribution; this can only be relevant either on very large scales or on very small highly non-linear 
scales |23j| - so ip = 4> to a good approximation for the range of wavelenghts of interest here. 


Collecting second order terms and adding to the first order term, we emerge with an expression for the perpendicular 
deflection accurate to second order in the gravitational potential: 


x±( 0 ) = -2 / dx{Xs-x) 


(V 1} (x {;o) (x, &j) + i</> (2) (£ (o) (x,0))) 
d 2 ^ 1 \x { -°\x,d)) [ x j (x^> (x'J)) 
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Note that the first set of terms on the right /2 is simply the full nonlinear potential out to second order; 

i.e., the metric (Eq. (1101) 1 contains this combination. So we might expect the full deflection to be sensitive to the 
fully nonlinear gravitational potential. Indeed this is the assumption built into all previous work. Our second order 
treatment justifies this assumption. Also note that the second term is indeed the Born correction alluded to in Eq. 0 
since twice the inner integral is equal to In fact, the so-called lens-lens correction is also included in this term: 
when differentiating x\ to get the shear, the derivative acting on <p" here gives what is usually called the Born 
correction, while the derivative acting on the <p' term inside the inner integral gives the lens-lens correction. We can 
rewrite the perpendicular deflection as 

(9) = ~2j dx{Xa-x)-^-^(x(xJ)) ( 22 ) 

where cf> and x now include all orders in perturbations: (f> = + (1/2 )(f>^ 2 ^ + ... and similarly for x. Eq. iPl) is the 

usual starting point for ray-tracing simulations. We have now rigorously justified the standard ray-tracing approach, 
at least to second order in geodesic corrections. 
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V. CONCLUSIONS 

We have jusitifed the standard ray tracing formula of gravitational lensing. This argument connects much of the 
formal work on second order perturbations in cosmology with the more phenomenological approach used in len sing 
studies. A similar result holds for the luminosity distance: the first order expression for the luminosity distance (24| 
can be used as the starting point for the second order result, without going back to the geodesic equation. 

Our primary result, Eq. (PI). holds only on sub-horizon scales. On very large scales, many of the corrections 
we were able to neglect are of the same order as the Born correction and the nonlinear evolution of the potential. 
However, on large scales, we expect all such corrections to be small, of order 10 -5 , since they are not enhanced by 
spatial variations. Therefore, our results are valid in the regime in which corrections are measurable. The one possible 
caveat to this claim is if super-horizon perturbations can influence the local observables (2^. 

This work is supported by the DOE and by NASA grant NAG5-10842. 
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